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1.  INTRODUCTION 


The  work  reported  here  was  done  at  the  Ballistic  Modeling  Diviaon  of  the  Ballistic  Research 
Laboratory,  (B^)  in  support  of  theoretical  and  experimental  efforts  througout  the  BRJL  The 
objective  of  this  work  was  the  refinement  and  maintenance  of  the  Eulerian,  hydrodynamics 
code,  BRLHELP^  at  the  state  of  the  an,  for  application  to  problems  in  both  basic  research  and 
experimental  studies.  This  objective  was  approached  through  the  analysis  of  calculations  of 
problems  of  current  interest,  and  of  the  applicability  of  various  models  in  BRLHELP  to  these 
problems. 

The  HELP^  code,  an  evolutionary  antecedem  of  BRLHELP,  has  been  used,  for  some  time, 
as  a  predictive  tool  in  warhead  mechanics,  and  intermediate  and  terminal  ballistics.  Its  use  is 
limited,  however,  by  serious  problems  with  the  accuracy  and  robustness  of  some  of  its  models. 
BRLHELP  is  devoid  of  most,  if  not  all,  of  these  problems. .  It  has  been  modified  to  calculate 
internal  energy  accurately^,  has  more  robust  treatments  of  certain  numerical  problems  inherent 
in  the  old  HELP  code,  and  has  temperature  available  as  an  output  variable  for  metals  and  alloys 
via  the  BRLGRAY^  equation  of  state. 


2.  BRIEF  DESCRIPTION  OF  MODELS 

Detailed  descriptions  of  the  mathematical  and  numerical  models  in  BRLHELP  can  be  found 
in  References  1-3.  A  brief  description  is  given  here. 

BRLHELP  is  a  two-dimensional,  Eulerian  hydrodynamics  code.  It  describes  unsteady 
multi-material  problems,  treats  material  strength  as  an  elastic-plastic  phenomenon,  and 
incorporates  a  multiphase  EOS  for  metals. 

The  code  uses  massless  tracer  particles  to  define  material  interfaces  and  free  surfaces.  These 
particles  are  placed  along  surface  boundaries  and  move  with  the  local  material  velocity.  This 
Lagrangian-like  description  provides  an  accurate  definition  of  the  material  interfaces  and  free 
surfaces  without  sacrificing  tlie  ease  of  treatment  of  extreme  distortion  by  the  Eulerian  model. 

Input  to  the  code  consists  of  an  initial  geometric  configuration,  initial  values,  explosive 
loading  parameters,  material  properties,  and  desired  editing  features. 

The  geometry  package  can  describe  arty  configuration  decomposable  into  straight  lines,  arcs 
of  circles,  and  arcs  of  ellipses.  The  high-explosive  detonation  model  can  treat  multiple 
explosives,  multiple  detonation  points,  and  the  presence  of  wave  shaping  devices.  Sliplires  may 
be  employed  at  material  interfaces  to  allow  different  materMs  to  develop  velocities 
commensurate  with  their  material  properties,  rather  than  a  common  interface  velocity. 


/.  Josej^i  Laceiera,  Janet  E.  Lacetera,  James  A.  Schmitts  T/ie  BRL  7600  Version  cf  the  HELP  CodJ*  Ballistic  Research  La¬ 
boratories  Report:  ARBRL-TR-02yi9,  Jan  SO.  (AD  # AO  82804) 

2.  Laura  J.  Hagpman,  et  al.,  '"HELP,  A  Multi-material  Eulerian  ProrjamJbr  Compressible  Fluid  and  Elastic-Plastic  Flows  in 
Two  Space  Dimensions  and  Timd",  Systems,  Science,  and  Software  Report  TR-76-45-BK2,  Apr  1976, 

i.  James  A.  SHimitt,  Tmneation  Error  Terms  in  the  Kinetic  Energy  Calculation  in  the  HELP  Algprithm  and  Their  Conse- 
0iencef  ,  J.  Cmnp.  Phys,  Vol  35,  No  2,  Apr  1990, 

4.  Joseph  Lacetera  '"BRLGRAY:  The  Ballistic  Research  Laboratory  Version  of  the  GRAV  Ectuotion  pfStat^,  ARBRL-TR- 
02253,  Aur  mo.  (AD  #A091246) 
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Output  includes  pressure,  density,  material  velocity  vectors,  specific  internal  energy,  and 
deviatoric  stress  components  as  ftinctions  of  space  and  time.  Temperature  of  metals  and  alloys 
is  also  available  when  BRLGRAY  is  chosen  as  the  equation  of  state. 


3.  APPLICATIONS 

BRLHELP  is  a  necessary  numerical  tool  for  the  perfonriance  of  optimization  and 
prefabrication  studies  for  materials  and  geometries  proposed  in  shaped-charge  warhead  desiga 
It  has  the  capability  of  modelling  the  collapse  of  the  liner  and  the  formation  and  developmem 
of  the  jet  and  the  slug.  It  provides  information  which  can  be  used  to  determine  jet-tip  velocity, 
collapse  velocity  and  collapse  angle,  liner  position  as  a  function  of  time,  stretching  rate  of  the 
jet,  and  strength  effects. 

The  HELP  code  had  been  used  with  success  in  predicting  jet-tip  velocities  and  liner  collapse 
for  a  variety  of  geometric  designs^”^.  However,  H^P  overestimated  the  internal  energy  in  the 
collapsed  portion  of  conical  shaped-charge  liner  materials.  This  overestimation  led,  in  some 
cases,  to  unreasonably  low  estimates  of  jet  densities^  Because  of  this  an  artificial  bound  on  the 
internal  energy  in  metals  in  shaped  charge  problems  was  installed  in  the  Tillotson*  EOS  in  an 
attempt  to  maintain  a  more  appropriate  jet  density.  Without  this  bound,  the  jet  tip  exhibited 
internal  energies  so  high  and  densities  so  low  as  to  imply  vaporizatioa 

In  BRLHELP,  these  problems  have  been  alleviated.  Internal  energies  are  calculated  quite 
accurately,  and  the  artificial  bound  on  internal  energy  in  metals  in  shaped-charge  liner  problems 
has  been  removed.  Calculations  are  performed  without  artificial  constructs  and  physically 
reasonable  results  are  being  obtained,  including  experimentally  verifiable  metal  temperatures. 

4.  UPDATED  VERSIONS  OF  THE  BRLHELP  CODE 

All  of  the  input  to  BRLHELP,  except  for  an  alphanumeric  header  card,  has  been  converted 
to  the  NAMELIST  format. 

The  NAMELIST  ii^jut  streams  are  called; 

OPTNS,  START,  RUN,  GRID,  MTRL,  GEOM,  SLIP,  and  DETN. 

In  the  HELP  code  input  variables  located  in  the  first  150  words  of  blank  common  are  called 
Z-variables  because  they  follow  Z(l),  the  first  word  of  blank  common,  and  therefore  can  be 
referenced  as  Z(I).  A  set  of  control  parameters  are  first  read  using  the  NAMELIST:  OPTNS. 
The  Z-variables  are  then  read  using  the  NAMELISTs;  START  aixi  RUN. 


5.  J.  T.  Harrison,  "Comparison  Between  the  Eiderian  Hydrodynamic  Computer  Code  (BRLSC)  and  Experimental  Collapse  of  a 
Shaped  Charge  LinetT.  BRL  Memorandum  Report  ARBRL-MR-02841,  Jun  1978  (AD  #059711), 

6.  Janet  E.  Lacetera  and  William  P.  Waiters.  "Theoretical  and  Experimental  Studies  of  Hemispherical  Shaped-Charge  Jjnerf, 
Transactions  of  the  25th  Conference  of  Army  Mathematicians,  ARO  Report  80-1,  Jan  1980. 

7.  A.  Klwan,  and  A.  Arbuckle,  "Study  of  Liner  Collapse,  Jet  Formation  and  Charaaeristicsfiom  Impulsive  Shaped-Charge  Sys- 
■  temf.  BRL  Report  No.  2028,  Nov  1977  (AD  #051342), 

8.  J.  H.  Tillotson,  "Metallic  Equation  of  State  for  Hypervdoctty  Impactf,  (Jeneral  Atomic  Report  No.  GS-3216,  Jid  1962. 
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The  fatal  error  caused  by  pure  cells  overemptying  has  been  corrected  using  coding 
generated  by  Janet  Lacetera’.  This  coding  is  discussed  in  Section  6  of  this  report  Pressure 
iteration  problems,  inherent  in  old  versions  of  the  code,  are  discussed  in  Section  7  and  a  new 
mass-wei^ted  average  pressure  scheme  is  discussed  in  ^tion  8.  A  call  to  a  subroutine  called 
BRLEOS  has  been  installed  in  the  old  subroutine  EQST,  in  place  of  equation-of-state  coding,  to 
allow  access  to  the  RRLBOS  code  system’®  \siiich  includes  BRLGRAY,  the  three-phase  EOS  for 
metals. 


There  have  been  several  versions  of  BRLHELP  active  on  the  BRL  7600.  They  are  cycles 
14  through  23  of  the  permanent  file  BRLHELP,  with  ID=CMCHYDRO.  They  differ  as 
follows: 


CY-14 

CY-15 

CY«=16 

CY=17 

CY=18 

CY=19 

CY=20 

CY-21 

CY=22 

CY-23 


BRL  7600  version  of  UNIV AC-1 108  HELP\ 

Cycle  14  modified  for  accurate  energy  calculation^. 

Cycle  14  with  complete  NAMELIST  ii^uL 
Cycle  15  with  complete  NAMELIST  ir^ut 

Cycle  17  with  OVERMT  correction  and  temporary  pressure  iteration  correctioiL 

C^le  18  with  FTN  compiler  modificatioa 

C^le  19  with  mass- weighted  average  pressure  for  mixed  cells. 

Cycle  20  with  most  of  EQST  functions  transferred  to  BRLEOS 
Cycle  21  with  remaining  EQST  functions  added  to  BRLEOS. 

C^le  22  with  temperature  map  for  metals. 


Cycles  14  through  22  have  been  retired. 


A  listing  of  the  jobsiream  which  created  cycle  17  of  BRLHELP,  after  cycle  16  had  been 
created,  is  shown  in  Reference  1.  Listings  showing  the  jobstreams  for  creating  later  cycles  are 
shown  in  Appendix  A  of  this  report. 


5.  DISCUSSION  OF  INPUT 

A.  CapaMities  and  Limitations 

The  problem  generator  allows  the  user  to  generate  material  interfaces  which  are  composed 
of  straight  line  segments,  or  portions  of  circles  or  ellipses,  with  the  followirig  limitations: 

1.  The  number  of  material  packages  in  a  calculation  is  limited  only  by  the  dimensions  of  the 
material  arrays. 

2.  A  single  cell  mn  contain  ar^  number  of  material  interfaces  and  ar^  number  of  distinct 
materials. 

3.  Different  material  packages  can  contain  the  same  material  with  different  initial  conditions 
(velocities,  energies,  densities). 

4.  A  package  can  be  divided  into  any  number  of  ^atially  distinct  subpackages. 


9.  Janet  Lacetera,  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground  Private  communication 

Id.  Jose/Jt  Lacetera  and  William  Buchheister,  "BRUVFF  and  the  BRLEOS  Code  Systenf,  Ballistic  Researdt  Laboratory  In¬ 
terim  Memoraruium  Report  No.  655,  Aug  1979. 
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5.  The  material  in  each  package  must  be  homogeneous. 

6.  The  grid  boundaries  do  not  act  as  free  surfaces.  The  code  assumes  that  a  package  that 
extends  to  the  edge  of  the  grid  is  infimte  in  that  direction 

7.  The  grid  must  have  at  least  three  rows  and/or  columns. 

8.  A  calculation  can  be  made  in  cylindrical  or  rectangular  coordinates. 

9.  The  top  and  right  grid  boundaries  are  transmittive.  The  left  boundary  is  reflective  (z-axis 
in  cylindrical  coordinates).  The  bottom  boundary  can  be  either  transmittive  or  reflective. 


B.  Edit  Control  Parameters  and  Z-Variables 

Edit  control  parameto's  are  read  using  the  NAMELIST:  OPTNS.  The  Z-variables  are  then 
read  using  the  NAMELISTS:  START  and  RUN.  These  are  defined  below. 

NAMELIST  OPTNS:  (Edit  control  parameters;  it  should  be  TX)ted  that 
the  first  three  of  these  parameters  suppress  diagnostic  print  and 
should  be  used  with  caution) 


IFLGST: 


IPMADJ: 


—0,  default 

— 1,  suppresses  message  "SUM  OF  FRACS  NOT  EQUAL  TO  TOTAL..." 
— 0,  default 

— 1,  suppresses  message  "EVACUATION  OF  MATERIAL  N..." 


ITPHSE: 


— 0,  default 

—  1,  suppresses  message  "MASS  EVAPORATED  DUE  TO  ROUNDOFF..." 


NPRTOP: 


“0,  defaidt 

=1,  suppresses  edit  of  sliplines 
-=2,  suppresses  long  edit  on  last  cycle 
—3,  combines  effects  of  (1 )  and  (2). 


TIMMAX:  The  maximum  time  in  CPU  secorrds  (decimal)  that  the  code  wU 

use  in  the  calculation 

For  NPRTOP  —2  or  3  TIMMAX  should  be  about  4  seconds  less 
than  the  jobcard  time;  otherwise  TIMMAX  should  be  about  10 
seconds  less  than  the  jobcard  time. 
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NAMELIST  START 


PK(1); 

Problem  identifier  (numeric) 

PK(2): 

Restart  cycle  number 

PK(3); 

-0 

for  startup  run; 

—1 

for  restart  without  short  edit  print  of  restan  cycle; 

—2 

for  restart  with  short  edit  print  of  restan  cycle. 

NAMELIST  RUN 

Contains  all  Z-variables  not  contained  in  NAMELIST  START. 


C.  Cell  Dimensions 

Ir^ut  data  which  define  cell  dimensions  are  read  into  NAMELIST  GRID.  Dimensions  of 
cells,  in  the  area  of  interest,  should  not  vary  by  more  than  10%,  nor  should  large  aspect  ratios 
(e.g.  dx/dy  >  2)  be  used  near  the  region  of  interest. 

D. 

Initial  Conditions 

In  order  for  the  material  tracer  particle  routine  ®  function,  the  flag  of  the  cell  in  the  bottom 
left  corner  of  the  grid,  MFLAG{2),  and  the  material  code  number  must  be  defined.  If  the  cell 
is  a  pure  cell  of  material  N,  set  MFLAG(2)— N;  if  it  will  be  an  interface,  then  set 
MFLAG(2)=-1;  if  it  will  be  a  void  cell,  set  MFLAG(2)“0.  (In  the  case  of  MFLAG(2)=-1  its 
value  ^  be  set  >  1(X)  after  the  interfaces  are  processed  by  CALFRC  and  VOLFND.) 

MFLAG{2),  material  code  numbers  (listed  in  EQST),  internal  energy  (ergs/gm),  and  radial 
and  axial  velocities  (cm/s)  are  all  read  into  NAMELIST  MTRL.  The  equation- of -state 
constants  are  defined  in  EQST  and  COMDIM.  Variations  to  the  initial  conditions  would  require 
modifications  to  FILGRD.  The  initial  density  RHOIN  may  be  different  than  the  material’s 
*  normal  density  RH021,  defined  in  COMDIM. 


E.  Strength  Constants 

Shear  yield  strength  and  tensile  strength  of  each  material  package  are  also  read  into 
NAMELIST  MTRL.  If  a  material  has  no  strength,  AMDM(N)  must  be  set  to  1;  the  shear 
yield  strength  constants  will  be  set  to  2ero  automatically  if  not  read  into  MTRL. 
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F.  Material  Tracer  Partides 


Material  package  boundaries  are  described  by  a  series  of  straight  lines  and/or  arcs  of  cirles, 
and/or  arcs  of  ellipses.  Ti^ut  is  read  by  TSETUP  into  NAMELIST  GEOM,  and  is  defined 
below.  Data  is  re^  into  GEOM  for  each  segment,  and  is  read  in  an  order  that  places  the 
package  to  the  left  of  the  segment.  Tracer  particles  are  not  placed  along  package  toundaries 
that  are  also  grid  boundaries. 


LT: 


—  1  straight  horizontal  line; 

—  2  straight  vertical  line; 

—  3  straight  diagonal  line; 

—  4  arc  of  circle 

-=  5  arc  of  ellipse 


(LT  must  be  negative  for  last  segment  of  a  package  or  subpackage; 
LT  must  equal  100  after  all  boundaries  have  been  defined.) 


MPN;  Material  package  number  (void  ■-  I'lMAT  +  1) 
NPTS:  Number  of  points  to  be  placed  along  line  or  arc 


XI S;  starting  X-coordinate 

YIS:  starting  Y-coordinate 

X2F;  final  X-coordinate 

Y2f;  final  Y-coordinate 

THl;  starting  angle  of  arc 

TH2;  terminal  angle  of  arc 

SLA;  X-coordinate  of  center  of  circle;  or  eccentricity  parameter  A 
in  equation  of  ellipse 

SLR  Y-coordinate  of  center  of  circle;  or  eccentricity  parameter  B 
in  equation  of  ellipse 
XCEN:  X,.  in  equation  of  ellipse 
YCEN;  Yc  in  equation  of  ellipse 
XR:  radius  of  circle 

Equation  of  ellipse: 

X-X.2  Y-Y/ 


An  arc  which  begins  parallel  to  the  positive  X-axis  and  ends  parallel  to  the  positive  Y-axis 
has  a  starting  angle  of  3.14159  and  a  terminal  angle  of  1.57080. 

LT  <  0  signals  TSETUP  to  generate  a  dummy  tracer  particle  \wth  coordinates  (-1000,0)  to 
mark  the  end  of  a  tracer  string  for  a  package  or  subpackage. 

The  segments  must  form  a  continuous  boundary  except  where  they  intersect  the  grid 
boundaries.  Note  also  that  each  interface  is  described  by  two  identical  sets  of  tracer  particles, 
so  the  number  of  tracers  specified  along  an  interface  must  be  the  same  for  each  package 
although  generated  in  opposite  directions. 

Generally,  there  should  be  at  least  two  tracers  per  cell  edge.  In  regions  of  great  distortion 
tracer  particles  may  become  sparse;  however,  activating  ADDTCR  (see  Z-variable  NADD) 
helps  to  ensure  that  the  tracer  particles  do  not  become  too  sparse.  The  2^variable  NTPMX 
limits  the  number  of  tracers  for  a  given  package  boundary.  NTPMX  should  equal  the  sum  of 
the  tracer  particle  arrays  TX  and  TY. 
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G.  Sipline  Endpoints 


If  sliplines  are  being  generated  (NOSLTP— 0)  SETUP  reads  data  defining  slipline  endpoirus 
into  NAMELIST  SLIP.  Each  material  is  specified  as  a  master,  slave,  or  neither.  All  of  the 
SLIP  variables  can  be  read  in  a  single  SLIP  stream.  These  variables  are  defined  below. 


MASTRD 

NSLAVD 

>IBGMD 

i^SD 

NENDMD 

NENDSD 


—  N  if  package  (of  material  N)  is  a  master 
“  0  if  package  is  a  slave 

“  N  if  package  is  a  slave 

—  0  if  package  is  a  master 

—  index  of  first  tracer  particle  of  the  slipline  for  master 

—  0  for  slave 

—  0  for  master 

•=  index  of  first  tracer  particle  of  the  slipline  for  slave 
•-  index  of  last  tracer  particle  of  the  slipline  for  master 
“  0  for  slave 
=  0  for  master 

“  index  of  last  tracer  particle  of  the  slipline  for  slave 


Sliplines  can  be  redefined  on  a  restart,  with  NOSLTP =0;  otherwise,  set  NOSLIP— 1  on  the 
restart. 


H.  High  Explosive  Detonation  Points 

SETUP  reads  data  defining  primary  atd  secoixJary  initiation  points  inm  NAMELIST  DETN. 
These  are  defined  below. 


IDET 


?roETS 

XDETS 

YDETS 

TIME 

XMIN 

YMIN 

XMAX 

YMAX 


—  1  for  primary  initiation  point 

—  2  for  secondary  initiation  point 
=  0  for  end  of  initiation  point  data 

—  material  package  numter  of  high  explosive 

—  X-coordinate  of  irudation  point 

—  Y-coordinate  of  initiation  point 

—  delay  time  if  any 

=  left  X-coordinate  of  search  area 
=  lower  Y-coordinate  of  search  area 
=  right  X-coordinate  of  search  area 

—  upper  Y-coordinate  of  search  area 


The  first  five  of  these  variables  are  read  into  DETN  once  for  each  initiation  point.  TIME 
corresponds  to  the  delay  time  if  any  for  primary  initiation  points.  The  secondary  points  must 
be  read  in  the  orto  in  w+iich  detonation  will  occur.  The  remaining  variables  are  read  inm 
DETN  once  for  each  detonation  point  to  defiine  a  search  area  for  the  code  to  determine 
detonation  times. 

A  card;  $DETN  IDET— 0$  must  be  read  in  problems  without  explosives. 
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I.  Restarts 


Restarts  require  only  the  NAMELISTs  OPTNS,  START,  and  RUN  unless  NOSLIP— 0  in 
which  case  NAMELIST  SLIP  is  also  read. 

On  restarts  NAMELIST  RUN  allows  the  user  to  redefine  ary  of  the  Z-block  •variables. 
However,  on  runs  which  terminated  on  TIMMAX  it  may  not  be  neccessary  to  charge  ary 
•variable  including  TSTOP  or  ICSTOP,  depending  on  whether  the  specified  •values  were  reached 
during  the  previous  rua 


6,  PURE  CELLS  OVEREMPTYING 

Subroutine  TPHASE  determines  whether  a  pure  cell  has  overemptied  by  checking  the  sign 
of  the  •variable  AMX(K),  where  AMX(K)  is  the  total  mass  in  the  Kth  cell  Previously,  a 
negati^ve  •value  of  AMX  would  trigger  an  error  exit.  The  current  codirig  conserves  mass  locally 
while  setting  the  mass  of  the  Kth  cell  to  zero.  Thus  the  code  is  allowed  to  continue  in  the  case 
where  the  overemptying  is  caused  by  small  numerical  error  during  a  stretching  condition 

In  anticipation  of  some  other  problem  causing  the  o^veremptying,  a  diagnostic  message  is 
printed  at  the  time  of  the  reco-very. 

The  modification  (IDENT  OVERMT)  is  sho'wn  as  part  of  the  UPDATE  ii^ut  stream  in 
Appendix  A. 


7.  PRESSURE  ITERATION  PROBLEMS  IN  OLD  VERSIONS  OF  HELP 

For  cells  containing  more  than  one  material,  the  old  HELP  code  determined  cell  pressure 
by  an  iteration  procedure.  The  densities  of  the  separate  materials  in  such  a  cell  were  -varied, 
with  specific  internal  energy  held  constant,  until  the  partial  pressures  con-verged  to  a  common 
-value.  The  procedure  was  di-vided  into  a  pre-iteration  calculation,  to  ensure  that  the  cell  was 
filled  exactly,  followed  by  the  iteration  calculation. 

For  the  mixed  cell  it  was  required  that  the  partial  pressures  satisfy 


Pi(Pi,Vi)-P(p,V) 


where  P(p,V)  is  the  cell  pressure. 

In  the  iteration  process,  the  Vj  were  -varied  until  the  Pi  -were  within  a  specified  tolerance  of  a 
common  pressure. 

Unfortunately,  the  pressure  iteration  routine  often  failed  to  converge  in  cases  where  the 
component  material  pressures  were  initially  approximately  equal,  but  the  relative  cell  -volumes 
were  quite  differeru. 
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One  reason  the  convergence  routine  failed  is  that  the  JWL”  EOS  has  a  pressure 
maxiitiuin  \sdiich  is  potentially  hazardous  to  such  a  pressure  iteration.  The  EOS  may  be  called 
during  the  iteration  process  with  a  pressure  above  this  maximum.  Such  an  occurrence  would 
cause  failure  to  converge.  To  avoid  this,  the  variable  RHOMAX  was  defined  in  a  DATA 
statement  in  subroutine  CDT  as  the  maximum  density  to  be  allowed  for  an  HE  during  the 
pressure  iteration.  (RHOMAX— 9  for  Comp  B.)  However,  this  procedure  itself  often  caused 
the  convergence  routine  to  fail. 

Temporary  codiitg  allowed  recovery  by  restoration  of  the  initial  values  of  masses  and 
pressures  in  the  component  materials.  The  values  for  the  material  which  occupied  the  largest 
volume  in  the  cell  were  then  taken  as  the  mixed  cell  values.  Consequently,  a  warning  message 
was  printed  indicating  pressure  iteration  failure  and  recovery.  Material  pressures  for  each 
material,  before  and  after  the  iteration  attempt,  were  also  printed.  This  coding  allowed  the 
study  of  mixed-cell  conditions  which  led  to  pressure  iteration  failure. 

Table  1  shows  output  from  an  actual  calculation  in  which  the  pressure-iteration  failure 
recoveiy  was  invoked.  Pressures,  before  and  after  iteration,  are  shown  for  three  consecutive 
calculational  cycles. '  The  pressures  for  each  material  are  approximately  equal  before  iteration, 
and  pressures  and  densities  for  one  or  both  materials  are  varied,  with  the  pressures  diverging 
during  the  iteration  process. 


TABLE  1 

Pressure  Iteration  Failure  Recovery  Data 


CY 

Pi 

Pi 

Pi 

Pz 

P\Ipi 

1055 

115.9 

2.25E-I-14 

8.66 

2.24E+14 

13.38 

121.9 

2.49E-1-14 

8.66 

2.24E-1-14 

1056 

115.8 

2.25E-1-14 

8.66 

2.24E+14 

13.37 

127.9 

2.75E-f-14 

9.56 

2.47E-I-14 

1057 

115.8 

2.25E+14 

8.66 

2.23E+14 

13.37 

133.8 

3.01E-I-14 

8.66 

2.23E+14 

11.  R  Lee,  M.  Finger,  andW.  Colbns.  'ML  Equation  qf  Stale 
Coeffcierasjbr  High  Explodve  d".  Lawrence  LivermoK  Laboratory. 
UCII>-161S9,  Jan  1973. 
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Table  2  shows  the  masses  and  mass  ratios,  at  5  cycle  intervals,  for  the  two  materials 
During  the  period  ^own,  material  1  is  passing  out  of  the  zone  and  material  2  is  occupying  an 
increasing  portion  of  the  total  cell  volume. 


TABLE  2 

Masses  and  Mass  Ratios  in  Mixed  Cells 


CY 

mj 

m2 

mi/m2 

1055 

6.67E-04 

8.59E4)2 

7.76&03 

1060 

5.42&04 

1.06E^1 

5.11E-03 

1065 

4.04E-04 

1.25E4)1 

3.23E-03 

1070 

2.68E-04 

1.42E-01 

1.89E-03 

1075 

1.10E04 

1.51E^1 

7.28E-04 

1080 

6.50E-05 

1.54&01 

4.22E4M 

1085 

3.34E-05 

1.55&01 

2.15E-04 

1090 

3.86E-06 

1.55E-01 

2.49E-05 

Qearly,  the  initial  pressure  of  materi^  2  at  the  beginning  of  the  process  is,  to  a  good 
approximation,  the  correct  mixed-cell  pressure.  First,  the  two  material  pressures  are 
approximately  equal  and  second,  material  2  occupies  a  dominant  ponion  of  the  cell  volume.  In 
such  a  case  the  iteration  process  itself  is  unnecessary  to  perform.  Furthermore,  it  is  apparent 
from  the  extremely  high  densities  of  material  1  that  the  iteration  construct  had  been  driving  the 
material  densities  to  urqjhysical  values,  well  before  the  failure  to  converge  had  indicated  a 
problem. 

Further  studies  wath  3-maierial  mixed  cells  led  m  the  conclusion  that  a  mass-weighted 
average  pressure  for  mixed  cells  was  the  more  efficient  approach.  The  coding  which  invokes 
the  new  average  pressure  calculation  is  discussed  in  the  next  section 
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8.  MASS-WEIGHTED  AVERAGE  PRESSURE  FOR  MIXED  CELLS 


Before  calculating  an  average  pressure  for  a  mixed  cell,  a  cell  filling  calculation  is  necessary 
because,  in  general,  the  masses  at  the  beginning  of  the  calculation,  do  not  exactly  fill  the  cell 
since  cell  masses  are  changed  in  the  TPHASE  portion  of  the  previous  calculational  cycle. 


The  code  computes  component  material  pressures  Pj  from  component  densities  p|  while 
holding  the  specific  internal  energies  Ej  constant  The  pi  are  then  varied  by  one  percent  to 
compute  new  P]  and  constant  energy  compressibilities: 


Ci^ 


APi* 


New  volumes  are  calculated  using 


AVi 


The  Pi  are  then  equated  to  l/Vj  where  the  total  cell  volume  is 


V-ZVi 


After  the  cell  filling,  the  mass-weighted  average  cell  pressures  are  computed  directly  from 


where 


and 


Ws  =  5>i 


where  the  mj  are  component  masses. 
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9.  FTN  COMPILER  ERROR 


The  currem  FTN  OPTl  compiler  at  the  HRL  computer  site  does  not  treat  certain 
NAMELIST  variable  addresses  correctly.  Apparently,  certain  three  letter  variable  addresses  are 
not  known  to  have  the  same  address  Iqr  the  NAMELIST  package  and  the  rest  of  the  code.  In 
particular,  values  of  the  parameters  RMU(I),  the  rigidity  moduli  of  the  materials,  the  UUR(l), 
initial  radial  (or  x)  velocities,  and  WA(T),  initial  axial  (or  y)  velocities,  read  as  non-zero  into 
the  location  reserved  for  them  by  NAMELIST,  remained  zero  to  the  rest  of  the  code  because 
the  variables  of  the  same  name  known  to  the  rest  of  the  code  had  different  addresses. 

For  the  duration  of  this  compiler  residing  on  the  RRL  Cyber  system,  corrective  coding 
(IDENT  MAY80FTNE),  shown  in  Appendix  A,  is  neccessary  for  avoiding  this  problem. 


10.  BRLEOS  CODE  SYSTEM 

In  old  versions  of  HELP  the  subroutine  EQST  performed  thermodynamic  calculations  using 
the  Tillotson  equation  of  state  for  metals  and  composites,  the  JWL  equation  of  state  for  explo¬ 
sives,  and  an  ideal  gas  equation  of  state  for  gases. 

In  the  new  version  of  the  code,  subroutine  EQST  is  used  to  initialize  certain  EOS  data  and 
call  subroutine  BRLEOS.  BRLEOS  includes  all  of  the  EOS  coding  of  the  old  EQST  plus  calls  to 
the  BRLGRAY  equation  of  state  and  the  BRLPUFF  equation  of  state. 

Applications  of  the  EOS’s  are  listed  in  Table  3  arxi  material  index  numbers  are  listed  in 
Table  4. 


TABLE  3 

Equation  of  State  Applications 

EOS 

Application 

Type 

BRLGRAY 

Metals/ Alloys 

Mixed-Phase 

Solid-Liquid-Vapor 

Tillotson 

Metals/Composites 

Joined  Compressed-Expanded 

Gamma-Law 

Ideal  Gas 

P={y-l]pE 

JWL 

Explosives 

P(E,p) 

BRLPUFF 

Metals/P  henolics 

Joined  Compressed-Expanded 
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TABLE 4 

Equation-of-State  Material  Irxlices 

Index 

Material 

EOS 

1 

Aluminum 

2 

Iron 

3 

Copper 

4 

Tantalum 

5 

Uranium 

6 

Magnesium 

7 

Berylium 

8 

Titanium 

9 

Tin 

BRLGRAY 

10 

Stainless 

11 

Niobium 

12 

Nickel 

13 

Molybdenum 

14 

Lead 

15 

Thorium 

16 

Tungsten 

17 

Gold 

18 

RHA 

101 

Tungsten 

102 

Copper 

103 

Iron 

104 

Aluminum 

105 

Berylium 

106 

Titanium 

107 

Nickel 

108 

Molybdenum 

109 

Thorium 

110 

Lead 

Tillotson 

111 

Chj 

112 

Granite 

113 

Andesite 

114 

Wet  Tuff 

115 

Dry  Tuff 

116 

Oil  Shale 

117 

Dolomite 

118 

Limestone 

119 

Halite 

120 

Ideal  Gas 

Gamma-Law 

121 

Comp  B 

122 

TNT 

123 

OctoI 

124 

PBX9404 

JWL 

125 

LX-14-0 

126 

Hivelite 

227 

Aluminum 

228 

Plexiglas 

BRLPUFF 

229 

Fused  Quartz 

230 

Quartz  Phenolic 

19 


The  EOS  constants  for  materials  112  to  119  were  given  by  Evans  and  Harlow^^  and  it  should  be 
noted  that  the  constants  in  the  Tillotson  EOS  for  the  materials  Granite  and  Andesite  are  valid 
only  for  pressures  less  than  150  kbar. 

11.  SAMPLE  CALCULATION 

A.  Sample  Input 

The  jobstream  corresponding  to  Sample  Problem  2  of  Reference  1  is  shown  in  Apperdix  R 
As  the  first  input  card  indicates,  this  calculation  is  of  a  copper,  hemispherical  liner  loaded  with 
unconfined  Comp  B. 

After  the  edit  control  parameters  and  Z-variables  are  read,  the  new  NAMELIST  iiput 
begins.  The  grid  that  is  specified  has  60  by  200  square  cells  of  .05cm  dimension.  The  next 
card  specifies  that  the  cell  in  the  bottom  left  corner  of  the  grid  is  a  pure  void  cell.  The  rest  of 
the  MTRL  input  specifies  material  1  as  copper  and  material  2  as  Comp  R 

The  GEOM  ii?)ut  begins  with  the  specification  of  an  arc  of  a  circle  (LT—4)  bounding 
material  1  with  300  tracers.  The  starting  angle  is  PT  and  the  terminal  angle  is  PT/2.  The  center 
of  the  circle  is  at  (0.0,6.35)  with  a  radius  of  1.905.  The  next  segment  is  a  horizontal  line  going 
from  (1.905,6.35)  to  (1.7907,6.35)  with  50  tracers.  The  description  of  the  liner  is  completed 
with  an  arc  (LT=-4)  of  radius  1.7907  corresponding  to  the  inner  radius  of  the  liner. 

The  remainder  of  the  GEOM  input  bounds  the  Comp  B  charge  (MPN=2)  and  the  void 
(MPN-3). 

In  this  problem  a  slipline  is  placed  along  the  interface  between  the  liner  and  the  explosive. 
The  SLIP  irqjut  is  straightforward.  The  indexing  begins  with  the  master  material  on  the  left 
starting  with  the  first  tracer  particle  of  the  master  material,  counting  up  to  the  last  master  tracer 
particle  on  the  interface  (300),  counting  back  along  the  interface  with  the  slave  material  on  the 
left  starting  with  the  first  slave  tracer  (301)  and  ending  with  the  last  slave  tracer  on  the 
interface  (600). 

For  this  explosive  configuration,  both  a  primary  and  secondary  detonation  point  are 
required.  They  are  placed  at  (0.0,2.54)  arxi  (1.82,5.62),  respectively.  Zero  delay  time  is 
specified  in  both  cases.  The  search  areas  specified,  for  the  code  to  calculate  cell  detonation 
times,  have  lower  left  and  upper  right  corners;  (0.0,2.54)  and  (2.159,6.35),  and  (0.0,5.59)  and 
(2.159,6.35),  respectively.  Less  overlap  would  be  desirable  but  the  computing  time  involved  is 
not  significant  What  is  important  is  the  placement  of  the  secondary  initiation  point  which 
allows  it  to  "shine"  on  the  ponions  of  the  explosive  not  covered  by  the  primary  point 


12.  M.  W.  Evans  and  F.H.  Harlow.  The  Particle-ln-Cen  Method Jbr 
Hydrodynamics  Calculation^',  Los  Alamos  Scientific  Laboratory 
Report  LA-2139.  Nov  19S7. 


20 


B.  Sample  Output 


Appendix  C  contains  the  output  for  Sample  Problem  2  run  on  BRLHEU,  cycle  23.  The 
first  page  shown  is  pan  of  the  dayfile.  A  copy  of  a  user-library,  LIBCMC,  containing  the 
BRU^S  codes  is  attached,  along  with  the  main  program  library  BRLHELP,  and  the 
MrXEDPHASEBOU?4DARY  data  tape  required  by  BRLGRAY.  This  was  a  shon  run,  using 
867.305  JOB  seconds  (corresponding  to  ’nMMAX-'900).  The  output  was  stored,  for  later  use 
on  a  restart,  on  a  permanent  disk  file  called  CUHEMT. 

The  first  page  of  actual  code  output  lists  the  edit  parameters.  The  riext  page  shows  the 
problem  identification  (read  on  the  header  card)  and  the  geometry  package  iirluding  definition 
of  the  slide-line  endpoints.  This  is  followed  by  a  description  of  the  explosive  package  on  two 
consecutive  pages.  The  Z-variables  are  printed  next,  followed  by  a  general  description  of  the 
problem  grid. 

The  material  in  package  1,  G  Copper,  is  copper  described  by  the  BRLGRAY  EOS.  The 
material  in  package  2,  J  Comp  B,  is  Comp  B  described  by  the  JWL  EOS.  Next  comes  the 
strength  constants  and  a  description  of  the  grid.  This  takes  several  pages  of  output  cotx;luding 
with  the  cpll-coordinates  of  the  tracer  panicles  for  each  material  package. 

The  compression  map  at  cycle  0  shows  the  compression  values  of  both  the  metal  liner  and 
the  explosive  (mixed-cells  are  indicated  by  asterisks)  to  have  the  symbolic  value  Z  Z  in  this 
case  has  the  value  of  1.0  where  in  general  the  value  of  Z  is  given  as  a  range  of  values,  for 
example  from  .980  to  1.026.  The  display  of  mixed  and  pure  cells,  shown  next,  shows  the  metal 
liner  as  material  1  and  the  explosive  as  material  2,  with  mixed  cells  denoted  symbolically  with 
the  letter  M.  Note  that  the  void  is  considered  as  a  material  for  computational  purposes. 

The  next  page  shown  is  the  end  of  the  cycle  0  printout  and  the  beginning  of  the  regular 
cyclic  printout  This  continues  until  it  is  time  for  a  tape  dump  and  printout  at  a  specified  time 
or  cycle,  or  at  the  end  of  the  calculatioa  In  this  case  the  next  tape  dump  was  at  the  specified 
time  of  2.0  pus  which  occurred  at  cycle  69. 

The  next  page  shown  is  the  pressure  map  at  2.0  /xs.  The  detonation  wave  has  not  quite 
reached  the  metal  liner,  and  the  BRLGRAY  EOS  has  not  yet  been  called  for  a  pressure  and 
temperature  calculatioa  When  the  shock  wave  does  hit  the  liner  and  the  EOS  is  called  for  the 
first  time,  it  prints  out  a  set  of  ii^ut  (ESTCON)  and  calculated  (EGONS)  equation-of-state 
constants.  This  occurs  during  the  regular  cyclic  printout 

The  last  page  of  output  shown  is  the  temperature  map.  Temperatures  in  the  liner  range 
symbolically  from  E  to  H,  indicating  temperatures  from  416  to  832  degrees  Kelvia  The  Q’s 
and  the  R  occurring  along  the  explosive/metal  interface  are  transient  values  which  have  little 
effect  on  the  calculatiort 
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12.  SUMMARY 


The  HELP  code  has  been  shown  to  be  an  excellent  tool  for  shaped-charge  liner  optimization 
studies;  however  there  have  been  some  nagging  problems  for  its  users.  Two  of  these  have 
been  addressed.  These  are  the  fatal  errors  caused  by.  (1)  over-emptying  of  pure  cells  and  (2) 
the  f^ure  of  the  pressure  iteration  routine  to  achieve  convergence  for  certain  mixed-cell 
situations.  A  cycle  23  version  of  BRLHELP  is  now  available  which  alleviates  these  problems. 
The  overemptying  of  pure  cells  has  been  converted  to  a  non-fatal  cortdition  and  pressure 
iteration  failure  for  mixed  cells  has  been  eliminated  by  employment  of  a  mass-weighted  average 
pressure  scheme. 

A  call  to  the  BRLEOS  code  system  has  been  installed  in  place  of  the  old  EQST  subroutine  to 
allow  access  to  a  variety  of  equations  of  state  including  BRLGRAY  which  provides  temperature 
as  an  output  variable.  A  problem  with  the  Fortran  compiler  on  the  BRL  Cyber  176  has  been 
circumvented. 

Ir^ut  to  the  BRLHELP  code  has  been  changed  to  the  NAMELIST  format.  Interpretation 
of  input  decks  is  thereby  less  cryptic  and  modifications  to  old  decks  less  prone  to  error. 

13.  RECOMMENDATIONS 

BRLHELP  should  be  considered  for  more  extensive  use,  at  the  BRL,  in  predlcti^'e  studies. 
Its  new  capability,  via  BRLGRAY,  to  produce  temperature  as  an  output  variable  should  prove 
most  useful  to  researchers  trying  to  determine  the  state  of  shocked  materials.  Its  new 
robusmess  should  prove  the  code  to  be  especially  useful  in  prefiring  analyses,  where  use  of  the 
code  can  save  considerable  expense  by  avoiding  trial  and  error  experimentation. 


JL  E.  Lee,  Finger,  andW,  Collins,  ^JWL  Equation  of  State 
Coefficients  for  High  Exj^osive  Lawrence  Livermore  Laboratory, 
UClD-16m,  Jan  1773, 

12.  M.  W.  Evans  and  F.H.  Harlow,  *The  Particle-in-Cell  Method  for 
Hydrodynamics  Calculationd* ,  Los  Alamos  Scientific  Laboratory 
Report  LA-2139,  ?Vov  1957. 
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APPENDIX  A;  UPDATE  LISTINGS 
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ICTRA,S1WFZ,T200.  UPDATE  TO  BRUJELP;  CY=17  TO  Cy=18 
ACCOUNT (MD440) 

REQUEST (B,  *PF) 

GETPF (A,BRLHELP,CY=17 , ID=CMCHYDRO) 

UPDATE (F,P=A,N=B) 

FIN  {I=COMPILE ,  ]>0 ,  PL=55000 ,  LCM=I ) 

CATALOG (B,BRLHELP,CY=8 , ID=CMCHYDRO) 

*EOR 

*IDENT  OVERMT 
*D  TPHASE.687,TPHASE.689 

C  ***  IF  PURE  CELL  OVEREMPTIED  CONSERVE  MASS  LOCALLY  *** 

C  ***  AND  ZERO  OUT  OVEREMPTIED  CELL.  *** 

1000  IF  {AMX(K))  1010,1070,1015 
1010  CONTINUE 

WRITE(6,1310)  I,J 
KKA=K-IMAX 

WRITE (6,1350)  I,J,K,MFIAG(K)  ,AMX(K)  ,AIX(K)  ,U(K)  ,V(K)  ,DETIM(K)  , 
1  MFIAG(KKA)  ,AMX  (KKA)  ,AIX(KKA)  ,U(KKA)  ,V(KKA)  ,DETIM(KKA)  ,ETH 
1350  FCRMAT(3I5,/I5,5E12.5,/I5,6E12.5) 

AMX  (KKA ) =AMX (KKA) -AMK (K ) 

AMX (K)=0. 

GO  TO  1060 
*IDENT  LACITERA 
*I  CDT.12 

DIMENSION  DNSAV(4) ,PRSAV(4) ,XMSAV(4) 

*I  CDT.71 

mSAV(L)=0. 

PRSAV(L)=0. 

XMSAV(L)=0. 

VMAX=0. 

PMXVL=0. 

*I  CDT.82 

E»3SAV(L)=RHO(L,M) 

XMSAV(L)=XMASS  (L,M) 

*I  CDT.129 

PRSAV(L)=PRS  (2,L) 

IF(VMAX.GT.VDL(L))  GO  TO  120 
VMAX=VOL (L) 

PMXVL=PRS (2,L) 

*D  CDT.246 

NK=260 

WRITE (6,560)  NK 
*D  CDT.249,CDT.250 
GO  TO  314 
*I  CDT.283 
301  CONTINUE 
*D  CDT.301 
NK=310 

WRITE (6,560)  NK 
*D  CDT.305,CDT.306 
314  DO  315  L=1,NMAT 
XMASS  (L,M)  =XMSAV(L) 
n^S(l,L)=r»^SAV(L) 

PRS(1,L)=PRSAV(L) 
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ICTRA,STMFZ,T200.  UPDATE  TO  BRIliELP;  CY'l?  TO  Cy-=18 
ACCOUNT  (MD440) 

REQUEST (B,  *PF) 

GETPF  {A,BRIiIELP,CY=17 ,  ID=CMCHYDRO) 

UPDATE  (F,P=A,N=B) 

FTO  {I»CCMPII£ ,  rr=0 ,  PI>55000 , 104=1 ) 

CATALOG  (B ,  BRIHEIP  ,Cy=8 ,  ID=CMCHYDRO) 

*EOR 

*IDENT  OVERMT 
*D  TPHASE.eeT^TPHASE.eSS 

C  ***  IF  PURE  CELL  OVEREMPTIED  CONSERVE  MASS  LOCALLY  *** 

C  ***  AND  ZERO  OUT  OVEREMPTIED  CELL.  *** 

1000  IF  (AMX(K))  1010,1070,1015 
1010  CONTINUE 

WRITE(6,1310)  I,J 
KKA=K-IMAX 

WRITE (6,1350)  I,J,K,MFIAG(K)  ,AMX(K)  ,AIX(K)  ,U(K)  ,V(K)  ,DETIM(K) 
1  MFIAG(KKA)  ,AMX(KKA)  ,AIX(KKA)  ,U{KKA)  ,V(KKA)  ,DETIM{KKA)  ,ETH 
1350  FORMAT (315, /I5, 5E12 . 5,/I5,6El2 . 5) 

AMX  (KKA )  =AMX  (KKA)  -AMK  (K ) 

AMX(K)=0. 

GO  TO  1060 
*IDENT  LACITERA 
*I  CDT.12 

DIMENSION  DNSAV(4) ,PRSAV(4) ,XMSAV(4) 

*I  CDT.71 

DNSAV(L)=0. 

PRSAV(L)=0. 

XMSAV(L)=0. 

VMAX=0. 

PMXVL=0. 

*I  CDT.82 

DNSAV(L)=RHO(L,M) 

XMSAV(L)=XMASS  (L,M) 

*I  CDT.129 

PRSAV(L)=PRS(2,L) 

IF(VMAX.GT.VOL(L))  GO  TO  120 
VMAX=VOL (L) 

PMXVL=PRS (2,L) 

*D  CDT.246 

NK=260 

WRITE (6,560)  NK 
*D  CDT.249,CDT.250 
GO  TO  314 
*I  CDT.283 
301  CONTINUE 
*D  CDT.301 
NK=310 

WRITE (6,560)  NK 
*D  CDT.305,CDT.306 
314  DO  315  L=1,NMAT 
XMASS(L,M)=XMSAV(L) 
ra^S(l,L)=I»ISAV(L) 

PRS(1,L)=PRSAV(L) 
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IF(PPS(1,L)  .GT.PMIN)  PSUM»1 
315  PAV=PMXVL 
*CCPy  CDT,CDT.303,CDT.304 
GO  TO  301 
*D  CDT.439 

560  FOPMAT(lHl,*PPESSURE  IlERATICN  FAILURE,  RECXJVERED  AT  CDT 

1  LINE  NUMBER*, 15) 


LCTRA,SIMFZ,T200.  UPDATE  TO  BRLHELP;  CY=18  TO  CY=19 
ACCOUNT (MD440) 

RECX3EST(B,  *PF) 

GETPF  (A ,  BRIiffiLP ,  ID=CMCHYDRO) 

UPDATE  (F,P=A,N=B) 

PIN  (I=COMPir£ , L=0 , PL=55000 ,  ICM=I ) 

CATALOG (B,BRLHELP, ID=CMCHyDRO,RP=999) 

*EOR 

*IDENT  MAYSOFTOE 
*D  GEOM.7 

1  CZERO,  STKl,  STK2,  STEZ,  TEMP,  AMDM 

*I  GECM.IO 

DO  100  IRMU=1,4 
100  RMU(IRMU)=TEMP(IRMU) 

C  ***  IHIS  CODING  NECESSARY  TO  AVOID  OPTl  COMPILER  ERROR  *** 


LCTRA,SIMFZ,T200.  UPDATE  TO  PRIRELP;  CY=19  TO  CY=20 
ACCOUNT (MD440) 

REQUEST (B,  *PF) 

GETPF  (A,  BRIRELP,  ID=CMCHYDRO) 

UPDATE  (F,P=A,N=B) 

FIN  (I=COMPir£ ,  L=0 ,  PL=55000 ,  ICM=I ) 

CATALOG  (B,BPIRELP,ID=CMCHYDRC,RP=999) 

*EOR 

*IDENT  CEIIRILL 
*YANK  lACITERA 
*D  CDT. 78 
*D  CDT. 95, CDT. 162 
C 

C  ***  COMPUTE  TWO  DENSITIES  FOR 
C  ***  EACH  MATERIAL  IN  CELL  USING 
C  ***  (1)  raiJSITIES  FRCM  LAST  CYCLE 

C  ***  (2)  THOSE  DENSITIES  INCREASED  BY  ONE  PERCENT 

C 

*OOPY  CDT,CDT.99,CDT.129 
*COPY  CDT,CDT.130,CDT.141 

IF  (VL(L).Gr.O)  GO  TO  170 
VL(L)=0.0 
XMASS (L,M)=0.0 
170  CONTINUE 


27 


*IDENT  PCALC 
*D  CDT.163,316 
*COPy  CDT,CDT.163,CDT.176 
*COPY  CDT,CDT.184,CDT.224 
C 

C  ***  MASS  WEIGHTED  AVERAGE  PRESSURE 
^  P(K)«=PSUM/WSUM 

C 

*COPY  CDT,CDT.308,CDT.316 
*IDENT  LCTSEP80 
*D  INPUT. 75, INPUT. 76 
G 


CALL  SETUP 

C  ***  CNLY  TO  READ  SLIP  LINE  DATA  AND  RETURN  *** 
C 

*B  SETUP. 13 
C 

C  ***  IF  RESTART  READ  SLIP  DATA  AND  RETURN  *** 
IF(PK(3).LT.O)  READ (5, SLIP) 

IF(PK(3) .LT.O)  RETURN 
C 

*D  GEOM.6,MAY80FTNE.l 
C 


NAMELIST/MTRL/  MFIAG,  MAT, 

1  CZERO,  STKl,  STK2, 

*D  MAY80FINE.2,MAY80PrNE.4 
DO  100  1=1,4 


RHOIN,  SSIEN,  UURX,  WAY, 
STEZ ,  RMUI ,  AMDM 


UUR(I)=UURX(I) 
WA(I)=WAY{I) 
100  RMU(I)=RMUI(I) 


C  ***  THIS  CODING  NECESSARY  TO  AVOID  OPTl  COMPILER  ERRQRf 

C  **  DETECTED  BY  JC^:  lACETERA  DURING  PERIOD  MAY-SEP  198( 
*I  SETUP. 11 


DIMENSION  UURX(4),  WAY (4),  RMUI  (4) 


ICTRA,STMFZ,T100.  UPDATE  TO  BRIHELP:  CY=20  TO  CY=21 
ACCOUNT  (MD440) 

REQUEST  (NEWPL,  *PF) 

GETPF  (OIDPL,ERIiffiLP,  ID=CMCHYDRD) 

UPDATE  (F,N) 

CATALOG  (NEWPL, BRIHELP,ID=CMCHYDRD) 

EXIT. 

*EOR 

*IDENT  EQSTEOS 
*D  BRL76.6 

1  TAPE10,TAPE1) 

*I  EOST.6 


*I 


COMMCN/EOSTEOS/XAMDM(4)  ,XEMIN,XGAMMA,XRDEPS,XMFK 
HELPCOM.71 

I£VEL  2,  PLWP,PLWC,PVOL,AIJPSV 
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*D  E0ST.11,E0ST.234 


XAMDM (MFK) «AMDM (MFK) 

XEWIN*EMIN 

XGAMMA^GAMMA 

XPOEPS-ROEPS 

XMFK=MFK 

ENERGY=AMAX1 (ENERGY ,0 . ) 

CALL  ERLEOS (ENERGY, RHOW,PRESUR,'rMPTR,CS, IV, N) 

*D  MEMEXPND.l 

COMMCN/MXCELL/  RHOZ(130) ,CNAUT(130) ,MAT(130) ,SAMMY(4,800) , 
*D  INPUT. 8 

DIMENSION  MNAME (130) 

*D  HELPEQ. 18, HELPED. 36 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


MATERIAL  PROPERTY  DATA  STATEMENTS 
MATERIALS  1  -  18  ARE  BRIGRAY  EOS  METALS 
MATERIALS  101-119  ARE  TILLCTSON  MATERIALS 
MATERIAL  120  IS  IDEAL  GAS 

NOTE  -  IDEAL  GAS  REFERENCE  DENSITY  RHOZ (120)  IS  THE 
NORMAL  DENSITY  OF  AIR.  THIS  VALUE  MUST  BE 
CHANGED  FCR  OTHER  GASES. 

EXPLOSIVE  121  IS  COMFB 
EXPLOSIVE  122  IS  TOT 
EXPLOSIVE  123  IS  OCTOL 
EXPLOSIVE  124  IS  FBX  9404 
EXPLOSIVE  125  IS  LX-14-0 
EXPLOSIVE  126  IS  HIVELITE 


DATA  (RHOZ(K) ,K=1,18)  / 


1 

2.70, 

7.87, 

8.94, 

16.65, 

19.08, 

2 

1.75, 

1.85, 

4.51, 

7.28, 

7.91, 

3 

8.58, 

8.86, 

10.22, 

11.36, 

11.68, 

4 

19.17, 

19.24, 

7.84/ 

DATA 

(RHOZ  (K) 

,K=101,119) 

/ 

1 

19.17, 

8.94, 

7.87, 

2.70, 

1.85, 

2 

4.51, 

8.86, 

10.22, 

11.68, 

11.36, 

3 

.90, 

2.70, 

2.70, 

1.97, 

1.70, 

4 

2.30, 

2.80, 

2.70, 

2.20/ 

DATA 

RHOZ(120)/.001293/ 

DATA 

(RHOZ (K) 

,K=121,126) 

/ 

1 

1.717, 

1.630, 

1.821, 

1.840, 

1.835 
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DATA  (CNAUT{K) 

1  4.01E5, 

2  4.78E5, 

3  2.89E5, 

4  3.49E5, 


K=101,119)/ 
3.95E5,  4.03E5, 

4.63E5,  5.15E5, 

2.58E5,  2.58E5, 

5.51E5,  3.85E5, 


5.27E5,  8.06E5, 
2.13E5,  2.03E5, 
2.24E5,  1.63E5, 
3.37E5/ 


*D  CDT.15,CErr.22 

GAMGAS(120)  =  GAMMA 

DATA  (GAMGAS(K),  K=121,126)/2.706,2.727,2.830,2.850,2.841,1.24/ 
DATA  (RHOMAX(K),  K=121,126)/9.0,  8.5,  9.5,  11.0,  9.5,  9.5/ 

C 


*D  DETIME. 7, DETIME. 18 
C 

C  EXPLOSIVE  121  IS  CCMPB 

C  EXPLOSIVE  122  IS  TOT 

C  EXPLOSIVE  123  IS  OCTOL 

C  EXPLOSIVE  124  IS  PBX  9404 

C  EXPLOSIVE  125  IS  LX-14-0 

C  EXPLOSIVE  126  IS  HIVELITE 

C 

C  EETVEL  IS  THE  DETONATION  VELOCITY 
C 


DATA  DETVEL(121)  /7.98E5/ 
DATA  DETVEL(122)  /6.93E5/ 
DATA  DETVEL(123)  /8.48E5/ 
DATA  DETVEL(124)  /8.80E5/ 
DATA  DETVEL(125)  /8.83E5/ 
DATA  DETVEL(126)  /.381E5/ 


*D  ADDENG.&,ADDENG.12 

DATA  (DENGY(K) ,K=121,126)  / 

2  4.95E+10,  3.68E+10,  5.27E+10,  5.54E+10,  1.02E+10,  4.20E+10/ 

C  EXPLOSIVE  121  IS  CCMPB 

C  EXPLOSIVE  122  IS  TOT 

C  EXPLOSIVE  123  IS  OCTOL 

C  EXPLOSIVE  124  IS  PBX  9404 

C  EXPLOSIVE  125  IS  LX-14-0 

C  EXPLOSIVE  126  IS  HIVELITE 

*D  INPUT.  10, INPUT.  16 

DATA  (MNAME(K) ,K=1,18)/ 

1  lOHG  ALUMINljM,10HG  IRON,10HG  COPPER, lOHG  TANTALUM, 

2  lOHG  URANIUM, lOHGMAGNESIUM, lOHG  BEPYLIUM,10HG  TITANIUM, 

3  lOHG  TIN,10HG  S-STEEL,10HG  NICBIUM,10HG  NICKEL, 

4  lOHG  MOLY,10HG  LEAD, lOHG  THORIUM, lOHG  TUNGSTEN, 

5  lOHG  GOLD, lOHG  RHA/ 


DATA  (MNAME(K),K=101,124)/ 
lOHT  TUNGSTEN, lOHT  COPPER, lOHT 

lOHT  BERYLIUM,10HT  TITANIUM, lOHT 


1 

2 

3 

4 

5 

6 
7 


lOHT  THORIUM, lOHT 
lOHT  POLYMERS, lOHT 
lOHT  DRY  TUFF, lOHT 
10HTLIMESTONE,10HT 
lOHJ  TNT,10HJ 


IRON, lOHT 
NICKEL, lOHT 

I£AD, 

GRANITE, lOHT  ANDESITE, lOHT 
OIL, lOHT  SHAIE,10HT 

HALITE, lOH  IDEAL  GAS,10HJ 
OCTOL, lOHJ  PBX,10HJ 


ALUMINUM, 

MOLY, 

WET  TUFF, 
DOLOMITE, 
COMP  B, 
9404/ 


30 


EftTA  MNAME(125)  /lOHJ  LX-14-0/ 

DATA  MNAME(126)  /lOHJ  HIVELIli:/ 

*D  INPUT. 99, INPUT. 107 

C  PRINT  INITIAL  CCNDITICNS 

WRITE(6,390) 

DO  60  L-1,NMAT 
MN=MAT(L) 

WRITE  (6,400)  L,RHOZ  (MN)  ,RHOIN (L)  ,SSIEN (L)  ,UUR  (L)  ,WA  (L)  ,MNAME  (MN) 
60  CONTINUE 
*D  ERL76.90 

400  PORMAT(I5,2F13.3,9X,1PE11.4,5X,1PE11.4,5X,1PE11.4,A10) 

*I  EDIT. 240 

305  CONTINUE 
*D  NOEDPRT.8 


LCTRA,SIMFZ,T200.  UPDATE  TO  ERIHELP;  CY=21  TO  CY=22 
ACCOUNT (MD440) 

RECUEST(B,  *PF) 

GETPF (TA , PLCMCEOS , ID=CMCeOS ) 

UPDATE (Q,P=TA,N=TEMPL) 

GETPF (A , BRLHELP , ID=CMCHYDRO) 

UPDATE (F, P=A  ,M=TEMPL,N=B) 

FTN  (I=COMPII£ ,  L=0 ,  PL=55000 ,  ICM=I ) 

CATALOG  {B,ERIiIELP,ID=CMCHYDRO,RP=999) 

*EOR 

♦COMPILE  BRIEOS 
*EOR 


*IDENT  ECSTEOSCR 
*D  ECSTEOS.45,ECSTEOS.47 

5  lOHT  DRY  TUFF, lour  OILSHAI£,10HT  DOLCMITE,10HTLIMESTONE, 

6  lOHT  HAUTE/ 

C 


C 


DATA  MNAME(120)  /lOH  IDEAL  GAS/ 


DATA  MNAME(121)  /lOHJ  COMP  B/ 

DATA  HNAME(122)  /lOHJ  ONT/ 

DATA  MNAME(123)  /lOHJ  OCTOL/ 

DATA  MNAME(124)  /lOHJ  PBX  9404/ 

♦IDENT  EOSEET 
*D  CDT.12 


DIMENSION  RHCMAX (130) ,  GAMGAS(130) 

*D  DETIME. 5 

1DETVEL(130) 

*D  ADDENG.4,ADDENG.6 

DIMENSION  DENGY(130) 

*D  ra:TIME.132 

IF(N.GT.120.AND.N.U.200)  GO  TO  105 
*D  EETIME.  135 

100  IF(MAT(M) .LE.120.QR.MAT(M) .GE.200)  GO  TO  120 
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*D  es:time.212 

IF(MAT(MFK)  .Gr.l20.AND.MAT(MFK)  .LT.200)  GO  TO  190 
*D  DETIME. 21 8 

IF  (MAT  (L)  .GT. 120.AND.MAT(L)  .LT. 200)  GO  TO  190 
*D  EETIME. 235 

IF(MAT(L)  .GT.120.AND.MAT{L)  .LT.200)  GO  TO  173 
*D  EETIME. 24 5 

IF  (MAT  (L)  .GT.  120. AND. MAT (L)  .LT.200)  GO  TO  176 
*D  DETIME. 256 

IF (MAT(MFL)  .Gr.l20.AND.MAT(MFL)  .LT.200)  GO  TO  181 
*D  DETIME. 261 

IF (MAT (L).Gr.l20.AND.MAT(L). LT.200)  GO  TO  181 
’►D  CETIME.272 

IF (MAT(MFL)  .GT.120.AND.MAT(MFL)  .LT.200)  GO  TO  185 
*D  EETIME. 276 

IF(MAT(L).Gr.l20.AND.MAT(L). LT.200)  GO  TO  185 
*D  IETIME.365 

IF(MAT(MFK)  .EE.120.OR.MAT(MFK)  .GE.200)  GO  TO  170 
*D  EETIME. 393 

IF(MAT(L)  .I£.120.OR.MAT(L)  .CE.200)  GO  TO  360 


LCTRA,S'IMFZ,T200.  UPDATE  TO  ERIHELP;  Cy=22  TO  Gf=23 
ACCOUNT (MD440) 

PECUEST(B,  *PF) 

GETPF  (A,BRIHELP,ID=CMCHYDPO) 

UPDATE  (F,P=^^,N=B) 

FTN  (I=COMPILE ,  L=0 ,  PL=55000 ,  ICM=I ) 

CATALOG  (B  ,BRLHELP,  ID=CMCHYDRO,RP=999) 

*EOR 

*IDENT  GEOMCR 
*D  GECM.15 

30  VJRITE (6,370)  XlS,yiS,X2F,y2F 
*IDENT  IMPOCTBO 
*I  HELPOOM.83 

COW1CN/TGRAY/TGRAY,  TXS  (12000),  TRS(4) 

LEVEL  2,  TGRAY,  TXG,  TRS 

*D  E0ST.3,E0ST.6 

C  ***  INITIALIZE  SCME  BRIGRAY  PARAMETERS  AND  CALL  BRIEOS 
C 

*B  ECeT.32 

TGRAY=0. 

*I  ECST.32 

TGRAY=TMPrR 
*B  CDT.23 

TGRAY=0 
*I  CDT.60 

TXG(K)=TGRAY 
*I  PCALC.12 

TRS (L)  =TGRAY 
*I  PCALC.35 
TSUM=0. 

*I  PCALC.54 

TSUM=AMAX1  (TSUM,TRS  (L) ) 
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ODD 


*I  PCALC.58 
C 

C  ***  MAXIMUM  lEMPEPATURE  FCR  ANY  METAL  IN  IWE  MIXED  CELL 
1XG(K)=TSUM 

*D  MAP.7 

DIMENSION  WSMAX  (6) ,  AI£(41),  XUM(41),  VALUE(41) 

*I  MAP. 71 

C  ***  TEMPERATURE 

DO  75  J=1,JDL 
DO  75  1=1, IDL 
K=(J-1)*IMAX+I+1 
M=IAES  (MFIAG(K)) 

N=MAT(M) 

IF(N.GE.IOO)  TXG(K)=0. 

IF(N.GE.100)  GO  TO  74 
IF(TX3(K)  .EO.O)  TXG(K)-300. 

74  CONTINUE 

75  WSMAX  (6)=AMAX1  (WSMAX (6)  ,ABS  (TXG  (K) ) ) 

PRINT  *,”IMAX=",WSMAX(6) 

*I  MAP. 139 

GO  TO  230 

***  TEMPERATURE 


225  iJ^JDL 

IF  (WSMAX  (NPROP)  .LE.O)  GO  TO  460 
WRITE (6,535)  CYCLE, T 

226  DO  228  1=1, IDL 
K=(J-1)*IMAX+I+1 

228  PROP(I)=T}G(K) 

*D  MAP. 234 

GO  TO  (90,120,150,180,210,226)  NPROP 
*D  MAP. 242 

GO  TO  (470,110,140,170,200,225,470)  NPROP 
*I  MAP. 259 

535  FORMAT (1H1,4X,15HTEMPERATORE  ,1X,6HCYCLE=,F8.1,4X,5HTIME=, 
*1PE11 . 5 , IX , 7HSECONDS//) 
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APPENDIX  B:  SAMPLE  INPUT 
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LCTIL,STMFZ,T1000,P0.  SAMPLE  PROBLEM  2:  CU  HEMISPHERE  LACETERA 
ACCOUNT (MD440) 

REQUEST  (TAPE7  ,  *P.F) 

GETPF (LIECMC,LIBCMC,  ID=CMCEOS) 

GETPF(A,  BRLHELP,  ID=CMCHYDRO) 

GETPF (TAPEl,  MIXEDPHASEBOUNDARY,  ID=CMCEOS) 

LIBRARY  (LIBCMC) 

EDITLIB. 


UPDATE  (F,P=A,N=B) 

FTN (I=COMPILE , L=0 , PL=55000 , LCM=I ) 
RFL  (160000) 

LGO. 

CATALOG (TAPE7,CUHEMI ,ID=CMCHYDRO) 
EXIT. 

*EOR 

LISTLIB (*, LIBCMC) 

*EOR 

*EOR 


SAMPLE  PROBLEM  2: COPPER  HEMISPHERE  LOADED  WITH  COMP  B  UNCONFINED 
$OPTNS  NPRTOP=2,  TIMMAX=900,  IFLGST=1 ,  IPMADJ=1,  ITPHSE=1 


$START 

PK(1)=1, 

PK  (2)=0 , 

PK  (3)=-0$ 

$RUN 

TSTOP=  6.0E-06 

,  PRDELT=  2.0E- 

-06, 

KUNITW=7 , 

KUNITR=7 , 

NOSLIP=0 , 

LVISC=1 , 

NMAT=2 , 

NMXCLS=800, 

NTRACR=5 , 

NTPMX=900,  REZ=0 

CYCPH3=-1 , 

NLINER=1 , 

NSLD=300 , 

MAPS=1 , 

IMAX=60 , 

JMAX=200 , 

DMIN=100 , 

IPR=100 , 
NFRELP=100, 

NDUMP7=1, 

IPCYCL=0 , 

11=3,  12=53, 

CVIS=-1 , 

IGM=0$ 

$GRID 

DX=60*.05,  DY= 

200*. 05$ 

$MTRL 

MFLAG  (2) =0 , 

$ 


MAT(1)=  3,  RHOINfl)=8.94, 

MAT(2)=  121,  RHOIN (2) =1 . 717 , 

CZERO(1)=2.35E+09,  STKl  (1 ) =6 . 95E+10 ,  STK2 (1 ) =5 . 50E+10 , 
STEZ(l)  =5.30E+09,  RMUI(l)  =4.55E+11,  AMDM(1)=  .9785, 

AMDM (2) =1 .0$ 
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$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GEOM 

$GECM 

$GEOM 

$GEOM 

$GEOM 

$SLIP 


$DETN 

$DETN 

$DETN 

$DETN 

$DETN 


LT=4,  MPN=1,  NPTS=300, 

TH1=3. 14159,  TH2=1. 57080, 

SLA=0.0,  ,  SLB=6.35,  XR=1.905$ 

LT=1,  MPN=1,  NPTS=50, 

X1S=1.905,  Y1S=6.35,  X2F=1.7907, 

LT=-4,  MPN=1,  NPTS=300, 

TH1=1. 57080,  TH2=3. 14159, 

SLA=0.0,  SLB=6.35,  XR=1.7907$ 

LT=1,  MPN=2,  NPTS=100, 

X1S=0.0,  Y1S=2.54,  X2F=2.159, 

LT=2,  MPN=2,  NPTS=200, 

X1S=2.159,  Y1S=2.54,  X2F=2.159, 

LT=1,  MPN=2,  NFTS=50, 

X1S=2.159,  Y1S=6.35,  X2F=1.905, 

LT=-4,  MPN=2,  NPTS=300, 

TH1=1. 57080,  TH2=3. 14159, 

SLA=0.0,  SLB=6.35,  XR=1.905$ 

LT=4,  MPN=3,  NPTS=300, 

TH1=3. 14159,  TH2=1. 57080, 

SLA=0.0,  SLB=6.35,  XR=1.7907$ 

LT=1,  MPN=3,  NPTS=50, 


Y2F=6 . 35$ 


Y2F=2 . 54$ 
Y2F=6 . 35$ 
Y2F=6. 35$ 


X1S=1.7907,  YIS 
LT=1 ,  MPN= 

X1S=1.905,  YIS 
LT=2 ,  MPN= 

X1S=2.159,  YIS 
LT=-1,  MPN= 

X1S=2.159,  YIS 
LT=100$ 
MASTRD(1)=1, 
NBGSD(l)  =0, 
MASTRD (2) =0 , 
NBGSD (2)=301  , 
IDET=1,  NDETS= 
IDET=2,  NDETS= 
IDET=0$ 


=6.35,  X2F=1 

3,  NPTS=50, 
=6.35,  X2F=2 

3,  NPTS=200, 

=6.35,  X2F=2 

3,  NPTS=100, 

=2.54,  X2F=0 

NSLAVD  (1)=0, 
NENDMD(1)=300, 
NSLAVD (2) =2 , 
NENDMD (2) =0 , 

2,  XDETS=0.0, 
2,  XDETS=1.82, 


.  905, 
.  159, 
.159, 
.0, 


Y2F=6. 35$ 
Y2F=6 . 35$ 
Y2F=2. 54$ 
Y2F=2 . 54$ 


NBGMD(1)=1, 

NENDSD(1)=0, 

NBGMD  (2) =0 , 

NENDSD (2) =600$ 
YDETS=2.54,  TIME=0.0$ 
YDETS=5.62,  TIME=0.0$ 


XMIN=0.0,  YMIN=2.54, 

XMIN=0.0,  YMIN=5.59, 


XMAX=2.159,  YMAX=6.35$ 
XMAX=2.159,  YMAX=6.35$ 
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APPENDIX  C:  SAMPLE  OUTPUT 
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I0/06/a0  SCOPC  2.1.5  B  R  L  VCR  0X1  ♦♦♦  11/06/80  80311 

STS  Devices  1l<>/  ^/PF  FLS-377K  FLL-1750K  f1XS»300K  MXL-1305K  MX8-1305B 


HH.Mf1.SS  CPU  SECOMD  aPICIM 
08.2''.28.MFA.  ♦♦ 
16.06.32  C0000.005  MFZ. 

16.06.32  00003.006  JOB. 

16.06.33  00030.026  JOB. 

16.08. 33  00000.032  JOB. 

16.06.33  00000.033  JIB. 

16.06.33  00003. 03A  JOB. 

16.06.33  00003.035  JOB. 

16.08.33  00000.03“'  MFZ. 

16.08.33  00000.003 
16.08.33  00000.00^ 
16.08.34  00000.02“^ 
16.03.34  00030.033 
16.08.34  00000.036 
16.08.34  00030.038 
16.08.34  00030.052 
16.00. 3S  00000.041  LHO. 

16.09.37  00000.059  USR. 

16.09.37  00000. 0-'S  USR. 

16.09,37  00000.073  USR. 

16.09.3“"  00000. 0''4  USR  . 

16,09,37  03000. 0'"4  LOD. 

16.09.37  00000.088  MFZ, 

16.09.36  00000.003 
16.09,36  00000.004 
16.09.4C  00030.027 
16.09.40  30033.033 
16.09.40  00000.036 
16,09.40  00000.038 
16.09,42  00000.087 
16.10.43  00334.098  USR, 

16.10.43  00004,009  L3D. 

16.17.47  00029.455  USR, 

16.17.47  00029.456  J3B. 

16.1"".50  00029, L10. 

16.25.46  00029.963  MFZ. 

16.25.48  00029.966  MFZ. 

16.25.48  00029.966  *JSR. 

16.26.00  00040.586  MFZ* 

16.26.00  00000,003  MFZ 
16.26.00  00030.004  JOB. 
16.26.02  00000.02'*  LOO. 
16,26.02  00000.033  '1FZ. 
16.26. J2  00030.036  MFZ. 


BP.L  rnS/BE  1.3  1499  VEP.  007 
H.CTIL#STMFZ#n000/P0. 

-ACCTUMTIMO***! 

“REQUESTCTAPE7/  ♦PF) 

-^ETPFCLXBCMCf  LIBCfIC/  IP-CMCEOS) 

-GET»FCA»  BP.LHEL®/  ID-C  MCHYDP.O) 

-GETPFITAPEl#  MIXEDPMaSE DCJUMDARY#  lO-CMCEOSJ 
-LIBRARY(LIBCMC) 

JM262  -  stage  PF  ST-MFZ  LFM-LIBCHC 


10/08/80 

SAMPLE  PROBLEM  2i  CU  HEMISPHERE  LACETERA 


MFZ 
JOB. 

LUD. 

MFZ. 

IFZ. 

MFZ. 

MFZ. 

-EOITLIB. 

F1RTRAM  LIBRARY  508 

snp 

042100  FINAL  EXECUTinN  FL. 

.014  CP  SECONDS  EXECUTION  TIME. 
-UPDATC(F/’*-A,M-B) 

JM262  ~  STAGE  PF  ST-MFZ  LFN-A 


-5JACP/STMFZ/T100/P7000.  ATT/CAT/PrG  SPOT  JOB 
-ACCTNITtCS***)  SYSTEM  JOB  5JACP 
-PFCOflV.  CAUL  ATTACM/CAT/PURGE  UTILITY 

PF646  -  PFMACP.n  -  ATTACH  -  LIBCMC  -  LIBCMC 
»F254  -  CYCLE  10  ATTACHED  FROM  StUSYSTEM 
JM262  -  STAGE  PF  SPOT  LFN-PF76C 

JM262  -  STAGE  PF  SPOT  LFf:.PF76C 


05/27/80 


MFZ, 
JOB. 
L  Tj. 
MFZ, 
MFZ. 
MFZ. 
MFZ. 


-SJACP/STMFZ#T100#P7000.  ATT/CAT/PRG  SPOT  JOB 
-ACC'TJJITICS***)  system  JOB  SJACP 
-PFCtTIV.  CALL  ATTACH/CAT/PUPCE  UTILITY 

PF646  -  PFMACRO  -  ATTACH  -  A  -  BRLHELP 

PF254  “  CYCLE  23  ATTACHED  FROM  SH-SYSTEM 
JM262  -  stage  PF  SPOT  LFN-Pf76C 

JM262  -  STAGE  PF  SPOT  LFN-PF76C 

MP3ATE  COMPLETED 

“Fine  I -COMPILE /L -0/PL- 55000/ LCM- I) 

25.34-"  CP  SECONDS  COMPILATION  TIME 
-RFL(160000J 
-LOO. 

LD610  -  FLS  REQUIRED  TO  LOAD  -  0043770  OU.COG 
L0603  -  EXECUTION  INITIATED  OS. EXP 
FORTRAN  LIBRARY  510  09/30/80 

J  J262  -  stage  PF  ST-MFZ  LFN-TAPEl 

-SJACP/ ST/1FZ/T100/P7000.  ATT/CAT/PRG  SPOT  JOB 
-ACCONMT(CS***I  SYSTEM  JOB  SJACP 
-PFCONV.  CALL  ATTACH/CAT/PUPCE  UTILITY 

PF646  -  PFMACRT  -  ATTACH  -  TA«»E1  -  MIXEDPHAS  EBQUfJO  ARY 
PF254  -  CYCLE  19  ATTACHED  FROM  SN-SYSTEM 


16.26.02 

00000 

.038  MFZ. 

JM262  -  STAGE  PF  SPOT 

IFN-PF76C 

16.26.02 

OOOJO 

.045  MFZ. 

JM262  -  STAGE  PF  SPOT 

LFfl-RF76C 

17.34.06 

C3097.272 

USP. 

EXIT 

17.34.06 

00897.272 

USR. 

:630'H3  FIMAL  EXECUTin'l  FL. 

17.34.06 

00897.2^3 

USR. 

967.305  CP  SECOfIDS  EXECUTION  TINE. 

17.34.06 

00397.276 

JOB. 

-CATALOG! TAP E7/CUIIEMI/  ID-CMCMYDROJ 

17.34.06 

00897.281 

MFZ. 

PF060 

-  CYCLE  3  CATALOGED  ON 

SN-SYSTEM 

17.34.06 

03897,282 

JOB. 

-EXIT. 

17.34.06 

00807,284 

MFZ, 

J  1166 

-  MAXIMUM  USER  SCM 

160000B 

WORDS 

17.34.  06 

0389*^.285 

MFZ. 

JM16T 

-  MAXIMUM  USER  LCM 

651000B 

WORDS 

17.34.00 

00897.285 

MFZ. 

JH70 

-  MAXIMUM  JS4I0  LCM 

2238 

BUFFERS 

17.34.06 

00897.205 

MFZ. 

RM'"’*0 

-  MAXIMUM  ACTIVE  FILES 

8 

17.34.06 

0089?, 265 

MFZ. 

R  1771 

-  3»EN/CLurE  CALLS 

48 

17.34.06 

00397. 285 

MFZ. 

RM77Z 

-  DATA  TRANSFER  CALLS 

193/876 

17.34.06 

00897. 2B5 

MFZ. 

RM7-"3 

-  contpol/positidnimg  calls 

1/002 

IT.  “JA..  n#. 

n'>o 07 . 

Mr  7 . 

-  ftM  HAT*  TDAurrrp 

1  n. 
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SAMPLE  PROBLEM  2  «  COPPER  HEMISPHERE  LOADED  WITH  COMP  B  ME  UMCONEIHED 
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DETOMATITN  time  CALCULATiaM  EXPLOSIVE  PACKAGES 
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PACKAGE  ‘HRMAL  INITIAL  CJNOITIONS 
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